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Abstract 

We study the probe of the gauge-phobic (or nearly gauge-phobic) heavy Higgs bosons (GPHB) 
at high energy hadron colliders including the 14 TeV LHC and the 50 TeV Super Proton-Proton 
Collider (SppC). We take the process pp —>■ titi, and study it at the hadron level including simulating 
the jet formation and top quark tagging (with jet substructure). We show that, for a CPHB with 
Mu < 800 CeV, can be determined by adjusting the value of in the theoretical pj,{bi) 
distribution to fit the observed Pt(&i) distribution, and the resonance peak can be seen at the SppC 
for Mh^SOO GeV and 1 TeV. 

PACS numbers: 14.80.Ec, 13.90.-|-i, 14.65.Ha 


I. Introduction 

The most important event in the 7 and 8 TeV runs 
of the LHC is the discovery of the 125 GeV Higgs boson 
[1]. The ATLAS and CMS collaborations have been mak¬ 
ing efforts to measure its couplings to the gauge bosons, 
T~^T~, and bb [3] [4] [5], and the obtained results are all 
consistent with the corresponding standard model (SM) 
couplings to the present experimental precision. How¬ 
ever, this does not imply that the SM is the final theory 
of fundamental interactions since the SM with a 125 GeV 
Higgs boson suffers from various shortcomings, such as 
the well-known theoretical problems [6] [7] [8]; the facts 
that it does not include the dark matter; it can nei¬ 
ther predict the mass of the Higgs boson nor predict the 
masses of all the fermions, etc. Searching for new physics 
beyond the SM is the most important goal of future par¬ 
ticle physics studies. So far there is no evidence of some 
well-known new physics models such as supersymmetry, 
large extra dimensions, etc. We know that most known 
new physics models contain more than one Higgs bosons 
in which the lightest one may be very close to the SM 
Higgs boson, and the masses of other heavy Higgs bosons 
are usually of the order of 10^ —10^ GeV. So that the dis¬ 
covered 125 GeV Higgs boson may be the lightest Higgs 
boson in certain new physics model, and probing other 
heavy Higgs bosons may be a feasible way of searching 
for new physics. Heavy Higgs bosons in the minimal 
supersymmetric extension of the SM (MSSM) and the 
two-Higgs-doulet model (2HDM) have been searched for 
at the LHC with negative results [9] . So that performing 
a model-independent search is more effective. 

In our previous paper [10], we proposed a sensitive 
way of probing anomalous heavy neutral Higgs bosons H 
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model-independently at the 14 TeV LHC via the process 
pp —7- VH* —> VVV —>• where V = W,Z. 

We showed that the resonance peak of H and the values 
of the anomalous coupling constants ^le 

measured experimentally provided the HVV couplings 
are not so small. Of course, this can not be applied 
to gauge-phobic (or nearly gauge-phobic) heavy Higgs 
bosons (GPHB) with vanishing (or very small) HVV 
couplings. For a GPHB H with mass in the 400 
GeV to a few TeV range, its main decay mode is 77 —>■ tt. 
If we simply consider the tt final state, the background 
will be extremely large. Ref. [11] showed that the process 
pp -T ttH tttt is a feasible process, and studied it at 
the parton level in the 2HDM. In this paper, we take the 
process pp —>■ tttt, and study the full tree level contri¬ 
bution at the hadron level including simulating the jet 
formation and top quark tagging (with jet substructure). 
We show that, with suitable kinematic cuts, the GPHB 
mass can be determined by adjusting the value of 
Mu in the theoretical Py(6i) distribution to fit the ex¬ 
perimentally measured py(6i) distribution for < 800 
GeV at both the LHC and the 50 TeV Super Proton- 
Proton Collider (SppC) considered in Beijing, and the 
resonance peak can be seen at the SppC for M^=800 
GeV and 1 TeV. 

II. Calculation and Results 

The general Yukawa coupling of the GPHB with the 
top quark can be written as 

Vti^Ht = C^yf^t^Ht, ( 1 ) 

where Ct is a parameter reflecting the deviation from the 
SM coupling. For simplicity, we take C^ k, In this 
paper, we take = 400 GeV, 600 GeV, 800 GeV, and 
1 TeV as examples to do the calculation. 

Ref. [12] shows that the next-to-leading-order correc¬ 
tion to the four-top production cross section in the SM 
is not so large. So we do the leading full tree level sim¬ 
ulation of pp -T titi at hadron colliders in this paper. 
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TABLE I: Detector acceptance according to DELPHES3 



fj- 

e 

jet 

photon 

l^lmarc 

2.4 

2.5 

5 

2.5 

PTmax (GsV) 

10 

10 

20 

0.5 


We use MadGraph5[13], FeynRules[14] and Pythia6.4[16] 
to simulate the signals and the backgrounds. We take 
CTEQ6.1[15] as the parton distribution function (PDF). 
DelphesS [17] and fastjet [18] are used to simulate de¬ 
tector acceptance and jet reconstruction. Our detector 
acceptance is shown in Table I referring to the design of 
CMS detector [19]. 

Some typical Feynman diagrams for the signal (S) and 
the irreducible background (IB) in pp —>■ tttt are depicted 
in Fig. 1. These two amplitudes will interfere with each 
other, so that they should be calculated together. 


or l£+2£“ {£ denotes p or e), and at least three tagged 
b jets. For reducible backgrounds (RB), there are three 
most important processes which can mimic the signal; 

• RBI: t~^t~bbti 

Processes with , bb^ and tt in the final state. One t 
(or t) decays semileptonically and the other i (or t) de¬ 
cays hadronically (cf. Fig. 2(a)). 

• RB2: ttt+jets and ttt+jets 

Processes with ti and another single top (anti-top) quark 
with extra jets in the final state. All the three top (anti¬ 
top) quarks decay semileptonically. We take one of the 
extra jets generated by matrix element and the other jets 
generated by parton shower (cf. Fig. 2(b)). 

• RB3: £+£“ + jets and bbf + jets 

Processes with £~^£~, bb, and a top (anti-top) quark with 
extra jets in the final state. The top (anti-top) quark 
decays semileptonically. We take one of the extra jets 
generated by matrix element and the other jets gener¬ 
ated by parton shower (cf. FIG. 2(c)). 



FIG. 1: Some typical Feynmann diagrams in pp —>■ tttt-. fa) 
for S, and (b) for IB. 

To suppress the SM background, we need some of the 
decay modes of the top quarks to include leptons. The 
CMS collaboration has studied the production of tttt in 
the SM with the final state including a single lepton and 
multiple jets at the 8 TeV LHC [20] , and has shown that 
it suffers from a very strong SM background. The AT¬ 
LAS Collaboration analyzed the production of tttt, and 
pointed out that, with the energy and integrated lumi¬ 
nosity of the LHC Run I, the most favorable final state is 
the one with two same-sign leptons [21]. Processes with 
fewer final state leptons can have larger cross sections 
but with lower signal to background ratios. In this pa¬ 
per, we shall present our simulation results for two kinds 
of final states. First we study the final states includ¬ 
ing three charged leptons (3-lepton mode) in which the 
SM reducible backgrounds (RBs) are highly suppressed. 
Next we study the final states including two opposite-sign 
leptons (2-lepton mode) in which the resonance peak of 
H can be observed. 

A. The case of 3-lepton mode 

We first study the case of 3-lepton mode, i.e. three 
of the four top (anti-top) quarks in the hnal state decay 
semileptonically while the other top (anti-top) quark de¬ 
cays hadronically. A signal event must contain 2£~^1£~ 



FIG. 2: Fenmann diagrams for typical RBs in the case of 
3-lepton mode. 

We generate events for S-I-IB, IB and RBs at the 14 
TeV LHC and the 50 TeV SppC, and do the detector 
simulation with DelphesS [17]. Anti-fer algorithm with 
i? = 0.5 are used to cluster jets and tag b jets. For 
the b-tagging efficiencies, we use the values in Ref. [22]. 
Then we apply the following kinematic cut: Requiring 
the event to contain three leptons (2£+lf“ or 2^“1^+) 
and at least three tagged b jets, i.e. 

A(£±) = 2, N(£^) = 1, N(b) > 3. (2) 

Since the S and the IB cross sections cannot be exactly 
separated, we define the IB, S, and total-background (B) 
cross section by = crg+iB(C't = 0), erg = crs+iB(C't 7^ 
0) — CTjb [23], and = CTjb -I- ct^b, respectively. Here 
is the total SM background. 

After this cut, the cross sections of S-I-IB, IB and RBs 
for various values of are shown in Table H. We see 
that RBs are negligibly small after the cut, and a smaller 





















3 


leads to larger signal cross section. Since the RBs 
are already negligible, we do not need to impose the top 


quark tagging requirement taking account of the jet sub¬ 
structure. 


TABLE II: Cross sections (in fb) of S-I-IB, IB and RBs for various values of after the cut in case of 3-lepton mode. 



400 GeV 

S-tIB 

600 GeV 

800 GeV 

1000 GeV 

IB 

RBI 

RB2 

RB3 

14 TeV LHC 

0.14 

0.10 

0.079 

0.066 

0.060 

1.3 X 10~“ 

1.3 X 10”'' 

6.4 X 10~' 

50 TeV SppC 

3.87 

3.31 

2.63 

2.21 

1.92 

0.039 

0.00032 

0.0039 


For an integrated luminosity Lint, the signal and back¬ 
ground event numbers are iVg = and 

respectively. For large enough and N^, the statisti¬ 
cal significance is defined as (Tgtnt = ^s/V^b- The 
needed for Icr, Scr, and 5cr deviations at the 14 TeV LHC 
and the 50 TeV SppC are shown in Table III. We see 


TABLE III: Integrated luminosity (in fb“^) needed for la, 3a 
and 5cr deviations at the 14 TeV LHC and the 50 TeV SppC 
for different values of M^. 



400 GeV 600 GeV 800 GeV 

1000 GeV 

la 

10 

37 

172 

1664 

14 TeV LHG 3cr 

87 

341 

1552 

14977 

5a 

241 

947 

4312 

41604 

la 

0.52 

1.0 

3.8 

23 

50 TeV SppG 3a 

4.6 

9.1 

34 

205 

5a 

13 

25 

96 

570 


AVg^B (bin)/N s_|_b — ANB(bin)/NB, with statistical er¬ 
rors for various values of with an integrated luminos¬ 
ity of 3 ab“^ at the 14 TeV LHC and the 50 TeV SppC, 
respectively. 




FIG. 3: AVs+B(bin)/Nsg_B — ANB(bin)/NB with statisti¬ 
cal errors at 14 TeV LHC: (a) Mfj—400 and 600 GeV, (b) 
Mh=800 and 1000 GeV. 


that the Scr case for = 600 GeV, and the 3a and 5 ct 
cases for Mj^ > 800 GeV at the LHG require very high 
luminosities which need the upgraded high luminosity 
LHC. All other cases can be reached at the present 14 
TeV LHC and the SppC. 

Since there are three missing neutrinos in the case of 
3-lepton mode, it is hard to construct the invariant mass 
distribution of the top quark pair from the decay of H 
to show the resonance peak of H. However, the values 
of may affect the distribution of certain kinematic 
observable from which we may determine the value of 
Mjj. In the hnal state, the b (b) quark is the secondary 
decay product which is more closely related to H than 
the charged leptons and ordinary jets (from W decays) 
do. Let bi be the b quark with largest prp. We choose 
the Pxibi) distribution to reflect the effects of different 
values of M^. Note that the unknown value of Ct also 
affects the Px{bi) distribution as an overall factor. To 
eliminate this effect, we separate the prp(bi) axis into 
a certain number of bins. Let AA(bin) be the num¬ 
ber of events within a bin at a certain value of prp(bi), 
and N = X)bin ^-^(bin) be the total number of events 
in the Py(6i) distribution. We then take the normal¬ 
ized distribution (ND), AV(bin)/N, for both S-l-B (in 
which the unknown Gj dependence is cancelled) and B 
(SM background). In Figs. 3 and 4, we plot the ND, 


E„ = 50TeV, l.|„, = 3 ab 
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FIG. 4: AVsg_B(bin)/Nsg_B — ANB(bin)/NB with statisti¬ 
cal errors at 50 TeV SppC: (a) M^=400 and 600 GeV, (b) 
Mh=S,00 and 1000 GeV. 


We see that the distributions are clearly distinguish¬ 
able for Mjj < 800 GeV at the LHC (the distributions 
for Mfj > 800 GeV can hardly be distinguished at the 
14 TeV LHC), and they are distinguishable for all val¬ 
ues of at the 50 TeV SppC. Thus the value of 
can be determined by adjusting the value of in the 
theoretical prp{bi) distribution to fit the experimentally 
measured prp{bi) distribution. 

Furthermore, if we adjust the value of in crg_|_B to 
fit the observed cross section, we may also determine the 
value of G( for the GPHB in nature. Of course the un¬ 
certainty depends on the experimental error. 
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B. The case of 2-lepton mode 



(c) 


FIG. 5: Feynmann diagrams for typical RBs in the case of 
2-lepton mode. 

To be able to show the H resonance peak, we consider 
the case of 2-lepton mode with the top quark pair from 
H decay decaying hadronically. When a top (anti-top) 
quark is highly boosted, its decay products can be clus¬ 
tered into a single fat jet by a clustering algorithm with 
a large R. And we can tag such a top (anti-top) quark 
at a high efficiency using jet substructure methods which 
are reviewed in paper [24] [25]. 

Some typical Feynman diagrams for S and IB before 
top quark decay are shown in Fig. 1. A signal event must 
contains two opposite-sign leptons, two tagged b jets, and 
two tagged t jets. The most important RBs which mimic 
the signal are: 

• RB-1: tijj 

Processes with a tt pair and two extra jets in the final 
state. The ti pair decay semileptonically (cf. Fig. 5(a)). 

• RB-2: £~^(.~bbjj 

Processes with a £'^i~ pair, a bb pair and two extra jets 
in the final state (cf. Fig 5(c)). In RB-1 and RB-2, 
we only generate two large-p-^ jets in the matrix element 
simulation in order to simplify the calculation. A fully 
merged background may modify the results, but the mod¬ 
ification will not be significant. 


• RB-3: £~^£~bbti 

Processes with a £'^£~ pair, a bb pair and a ti pair in the 
final state. The ti pair decay to jets (cf. Fig. 5(b)). 

The simulation is similar to the case of the 3-lepton 
mode, while, in addition, we should take the top quark 
tagging considering the jet substructure. We make jet 
clustering using the Cambridge-Aachen (CA) algorithm 
[26] with radius R = 1.2. The jet pruning algorithm [27] 
with parameters Zcut=0.I and RFactorcut=0-5 is applied 
to the CA jets for further suppressing RBs. If the trans¬ 
verse momentum of a pruned CA jet is larger than 350 
GeV, we retain it as a boosted jet. The components of 
the CA jets with pj, < 350 GeV are then clustered again 
by the anti-fer algorithm with R = 0.5. A b-tagging 
scheme with the same parameters as in Ref. [22] is also 
performed. We then impose the following cuts: 

Cut 1: Requiring the event to contain two isolated 
opposite-sign leptons and two tagged b jets, i.e. 

N{£+) = 1, N{£-) = 1, N{b) = 2. (3) 

Cut 2: Let be the jet with largest pj. and J 2 be the 
jet with the second largest prp. We require 

p^(jJ> 350GeV, |r?(jJ|<2 (4) 

This can pick up the events containing two boosted jets. 

Cut 3: Requiring the mass of the two boosted jets 
and fy to be in the neighborhood of m^, i.e. 

150 GeV < M(jfy < 220 GeV (5) 

This makes ji and fy two tagged boosted top jets. 

Cut 4: Since H is produced by top-quark fusion, its mo¬ 
mentum should not be so large. When the two boosted 
jets are the decay products of an s-channel H, the abso¬ 
lute value of the vector sum of their 3-momenta should 
be relatively small. So we require 

IpO'i)+p 02)I < ITeV. (6) 

This makes the two boosted top jets coming from H de¬ 
cay. 

The cross sections (in fb) of S-I-IB, IB, and three RBs at 
the 50 TeV SppC after the cuts are listed in Table IV. We 
see that after the four cuts, the backgrounds are reduced 
to the same order of magnitude as S-I-IB. 


TABLE IV: Cross sections (in fb) of S+IB, IB and RBs for various values of Mfj at the 50 TeV SppC after each cut in the 
case of 2-lepton mode. 



400 GeV 

S+IB 

600 GeV 

800 GeV 

1 TeV 

IB 

RB-1 

RB-2 

RB-3 

Cut 1 & 2 

1.05 

0.70 

0.65 

0.62 

0.40 

34.6 

1.79 

0.06 

Cut 3 

0.086 

0.061 

0.063 

0.07 

0.036 

0.30 

0.0096 

0.00049 

Cut 4 

0.032 

0.022 

0.022 

0.026 

0.013 

0.067 

0.0012 

0.00011 
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Since the two top jets from H decay for < 800 
GeV can hardly satisfy cut 2, we take the cases of = 
1 TeV and = 800 GeV to plot the invariant mass 
distributions for an integrated luminosity of 3 ab“^ at 
the 50 TeV SppC in Fig. 6. the resonance peaks can be 
observed. Note that the resonance peaks in Fig 6 are not 



luminosity of 3 ab“^ at the 50 TeV SppC: (a) Mfj = 1 TeV, 
(b) Mff = 800 GeV. The SM background (blue-dotted) is also 
plotted for comparison. 

clear enough for precision measurement of A more 
sophisticated, while complicated, top tagging algorithm 
[28] [29] may improve the results. 

III. Summary 

In this paper, we have studied the probe of the gauge- 
phobic heavy Higgs bosons H, usually appear in new 
physics models, at the 14 TeV LHG and the 50 TeV SppG 
via the process pp —?> titi mainly contributed from gluon 
fusion. We take the general (model-independent) Htt 
Yukawa interaction form in Eq. (1). 

Our calculation is at the hadron level including sim¬ 
ulating the jet formation and top quark tagging (with 
jet substructure) taking account of the requirement of 


the detector acceptance. For suppressing the SM back¬ 
ground, we take two of the decay modes of the top quarks 
including leptons, namely the 3-lepton mode and the 2- 
lepton mode, the former can have negligible SM RBs, 
and the latter can make the resonance peak of H visible. 

In the 3-lepton mode, the needed integrated luminosi¬ 
ties for Icr, 3cr and 5a deviations from the SM background 
are shown in Table IIL The case of 5cr for = 600 
GeV, and the cases of 3cr and 5cr for M^j > 800 GeV at 
the LHC require very high luminosities which need the 
upgraded high luminosity LHC. All other cases can be 
reached at the present 14 TeV LHC and the SppC. We 
have further shown that the mass of H can be experimen¬ 
tally determined by measuring the p,p(6i)-distribution at 
the 50 TeV SppC for an integrated luminosity of 3 ab“^ 
(cf. Figs. 3 and 4). 

In the 2-lepton mode, we imposed a series of cuts to 
suppress the background, and extracted the process that 
the t and t with hadronic decays are from the H decay. 
The plotted invariant mass distributions AI(j]^, at the 
50 TeV SppC for an integrated luminosity of 3 ab“^ are 
shown in Fig. 6 which shows that the resonance peak of 
H can be observed. But the peak is not clear enough to 
give a very accurate measurement for M^j. More sophis¬ 
ticated top tagging algorithm [28] [29] compared with the 
simple one used here may help to highlight the resonance 
peaks. 
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